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It is well known that the convection of liquid solution or melt liquid is inhibited by the magnetic 
field experimentally, and so that we can make the larger and better crystals at the cost of slow growth 
rate [lj, 0]. We shall show here another two effects due to the magnetic field. First, when the 
ionised molecule diffuses on the surface of crystal under the magnetic field applied normal to the 
surface, its diffusion constant decreases due to the magnetic effect. Second, even if the molecule 
has no electric charge, if it has electric dipole moment, the kinetic energy decreases by different 
mechanism. 
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I. CHARGED PARTICLE DIFFUSION 

Let us consider the motion of charged particle on flat 
surface under the magnetic field. Without magnetic field, 
there are various driving forces to make a motion of par- 
ticle, such as, random force due to temperature, fluid 
mechanical force, e.t.c. We introduce the velocity due 
to these forces by v. 



v = jj,F e ., 



(1) 



where u is the mobility. Then the electric current is 
given by 

J tot = env + J mag - (2) 

Where e is the electric charge, n is the particle density, 
and Jmag is the current induced by magnetic field, such 
as, 



J,, 



(—)(Jtot x B), 

en 



(3) 



where a is electric conductivity. By mixing these two 
equations, we obtain 

J tot = env + { — ){Jtot x B). (4) 
en 

From the iteration method, we obtain 



Jtot = env + ( — ){J to t x B) 
en 
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= env + a(vx B) - ( — )B 2 v + 0{B 3 ). (5) 
en 
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The third term means the inhibition of motion due to 
the magnetic field. 

To solve the full equation without approximation, we 
set the coordinates as, 



Jtot = (Jx,J y ,0), B = (0,0,B). (6) 
Then we have 



1 


-oB ' 




J X 




v x 


oB 


1 




Jy _ 


= en 


. v v . 



(7) 



where a — a /(en). By multiplying inverse matrix in 
both hand sides, we obtain 



Jx 



1 + a 2 B< 



(v x +aBv y ), J x 



l + a 2 B 



v y -aBv x ). 
(8) 



If we use mobility instead of electric conductivity, 



a = e nfi, 



(9) 



we obtain 

Jtot 



1 + (efiB) 2 



, n V X B . . 

^ITW' (10) 



Next we consider the case of diffusion. We simply ex- 
change the variable as 



D\7n = nv. 



(11) 



Then we have 



eD 



1 + (euB) 2 

eD 
1 + (euB) 2 



(d x n + e[iBd y n), (12) 
(d y n — e[iBd x n). (13) 
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We obtain from ([20 



X 2 o_ g 



FIG. 1: Electric dipole moment 



By using the charge conservation law 



- = div(Jtot), 
at 

we obtain the diffusion equation. 



(14) 



leu - 
ft = —R x B 
ml 



ml 



(22) 



By multiplying ft in both hand sides of (|22p . we obtain 



(23) 



2eB ■ 

9 = -(Ysin9 + Xcos9), 

ml 

where 

ft = (cos 9, sin 9), ft = (— sin 9, cos 0)0, 
n = (-sin#, cos 9)9 > — (cos6>, sin#)# 2 , 

and -R = (X, F) are utilised. 



r V 2 n 



1 + (e^B) 2 ' '" dt 
This shows the new effective diffusion constant. 

D 

D -> £»' = 



1 + (e/iS) 2 



(15) 



(16) 



In this model, energy dissipation is induced by the ex- 
istence of electric conductivity or mobility. The energy 
loss is easily calculated as, 

^-^(en)f-i,f } ^-^^E. (17) 



II. MOTION OF A PARTICLE WITH 
ELECTRIC DIPOLE MOMENT UNDER 
MAGNETIC FIELD 

Next we consider the motion of electric dipole moment 
which is planar to the surface under the magnetic field ap- 
plied normal to the surface. The distance of two charges 
are taken to be constant length /. Then the relative and 
center of mass coordinates are 



r = x± — X2 — In, R 

respectively. 

The motion of equations are 



X\ + x 2 



(18) 



mx\ = exi x B + nf, mx2 = ~ex2 X B — nf, (19) 



where / is the constraint force. 

From these two equations, we obtain 

•j. 2e u - 2/ _ 
f = —R x B + —n, 
m m 

R = —fx B. 

2m 

First we consider the relative motion. 



(20) 
(21) 



On the other hand, by multiplying ft in both hand sides 
of ([22]), we obtain 

Jill . 

/ = -{— 9 2 +eB{Ycos9-Xsin9)}. (24) 

This equation simply determines constraint force /. 

Next we consider the motion of center of mass. From 
(l2"Tj) . we obtain 



O eBl A 

X = 9 cos ( 

2m 

The time integration gives 
eBl 

X = — sin9 + U x , Y 
2m 



eBl ■ 
Y = —9sm9. 
2m 



(25) 



eBl 
2m 



cos9 + U y , (26) 



where U x , U y are the integration constants. 
By putting (f2T)|) into (|2^)l . we obtain 
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2eB 



ml 

Now let us denote 



■(Uy sin 9 + U x cost 



U x = U cos(j>, U y = Usmcj). 
Then we have equation 



2eB, 
ml 



-Ucos( 



(27) 



(28) 



(29) 



This is nothing but the equation of pendulum. 
To see it clearly, we write down the effective energy of 
this local system by using p = el. 



E eff = — 9 2 +pBU{sm(9 



+ 1}. (30) 



Compare to the pendulum weight, the lowest point cor- 
responds to the stable fixed point, which is shown in fig- 
ure 2 (left). On the other hand, the highest point cor- 
responds to the unstable fixed point. Note that when 
B < this relation is exchanged. 
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FIG. 2: Fixed points 

III. ENERGY CONSERVATION AND 
MECHANISM OF DISSIPATION 

Let us calculate the kinetic energy of center of mass. 



Ekin — 



2 V 



Y 



f ((P^l sine + U x f + (- ] § cos9 + U y f) 



{ P Bf 
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pB(s'me U x — cos 6* Uy) 



where M 



2M ' 2 ' * K ~* v ' 
const. + pBU sin(e - 0), (31) 

= 2m, the mass of dipole moment. Therefore 
dE k 



By using 



dt 

we obtain 



pBUe cos( 



dt 



p 2 m d e 2 



We define the moment of inertia / as 



Then we have 



dt (E km 



Ml 2 



= 0. 



(32) 



(33) 



(34) 



(35) 



The total energy is the sum of kinetic energy and rotat- 
ing energy that conserves. The key point here is, when 
the magnetic field is absent, each energy conserves in- 
dependently. But when the field applied, relative and 
center of mass motion couples. Therefore when the rela- 
tive motion is put in heat bath and give its energy to the 
bath, total energy decrease and we find the dissipation 
not only the relative motion but also the motion of center 
of mass. This mechanism gives the effect of dissipation 
to the motion of center of mass only when magnetic field 
applied. To realise the idea, we give a model in the next 
section. 



IV. MODEL OF DISSIPATION 

We introduce the dissipation only to the relative mo- 
tion, and see how the kinetic energy of center of mass 
decreases when the magnetic field applied. 



Then our equation for the relative (rotating) motion is 
described by 



— R x B 

m 



2/^ a u 
— n — -r 
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instead of (|2"0|) with positive a. 
And that gives 



2eB 
ml 



U cos( 



a6, 



instead of 

Then we obtain the interesting relation. 



dE, 
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;(E ki . 



A2 ) = -aie 2 < 0. 



(36) 



(37) 



(38) 



dt dV^' 1 ' 2' 
That is, E e ff is the dissipating total energy. From this 
dissipating property and from the form of E e ff, we see 
e — > (j) =p 7r/2 when t — > oo. The signature corresponds 
to the one of magnetic field, i.e. B = ± | B \. Then the 
kinetic energy E^m takes minimum value at t — > oo. 
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(39) 



Therefore the kinetic energy decreases due to the dissipa- 
tion of relative (rotational) motion only when the mag- 
netic field applied. Furthermore, it should be noted that 
the minimum kinetic energy is Non zero unless 

U = ujJ, (40) 

holds. Where the Lamor frequency lo c is defined by 

elBl (41) 



M 

From the existence of specific length / we have a constant 
lo c I with dimension of velocity . Then some kind of 
"resonance" occurs with velocity U that reduces the 
minimum kinetic energy to be zero. 

The simulation of the motion of electric dipole is given 
in 3 figures, (figure 3-5) Parameters are taken as U = 3, 
4> = 7T, e(t = 0) = 6{t = 0) = 0, uj c I = 7T, and a = 
0.4. Figure 3 shows the time dependence of (effective) 
total energy. Figure 4 shows the time dependence of the 
kinetic energy with its mean value for each one period, 
and Figure 5 shows the trajectory of the center of mass. 

From these result, our scenario that the kinetic (mean) 
energy dissipates when the magnetic field applied, is re- 
alized. 



V. THERMAL EFFECT TO OUR DIPOLE 
SYSTEM 

Next we put the system of relative motion in the heat 
bath. Then since the energy transfers from the rotat- 
ing motion to the heat bath, we may expect the similar 
dissipative dynamics occurs in the diffusion system. 
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FIG. 3: Total energy 




FIG. 4: Kinetic energy of center of mass 



By taking the ensemble average of equation (|26|). we 
have 



< v x > = 



P\B 
M 
p\B 

<%> = m~ 



<sm6>+U x , (42) 
< cos 9 > +U y . (43) 



This average is re-interpreted as the thermal average 
due to the Ergode hypothesis as follows. 



< f(0) >-- 



Jd9j dP f{9) exp[-(3H(9,P e )} 



Jd9j dP 9 exp[-PH(9,P e ) 
with Hamiltonian 



(44) 



for B > 0. And 



< sin 9 >~ cos ( 



< cos# >: 



(47) 



for B < 0. 
Therefore 

<v at > = -^t?^ cos <f> + U x = {1 - £JAi}Ekj(48) 
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The mean velocity becomes zero when 



P\B\ 
MU 



,1 



U 



= 1, 



this is related to the case Ekm = as shown in equation 
(f3U)) . This is NOT a dissipation, but is some "resonance" 
as we discussed already. 

The finite temperature case can be obtained by explicit 
calculation as well. 



< sin 9 > 



J de sin(— J + e + <j>) exp[— a(l — cos e)] 
J deexp[— a(l — cose)] 
J de cos e exp [a cos e] 



cos 



J dsexp[a cos e] 

d f 
cos 4>— log{ / cfeexp[a cose]}. 



(50) 



with 



a = /3BpU. 



In the same way, 



H = -^-Pl + pBU {sin{ 
mU 



+ !}• (45) 



In the zero temperature limit (3 — > oo, we have the 
lowest potential energy, that is, 



< sin 9 >~ — cos < 



< cos 9 >~ sin</>, (46) 
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FIG. 5: Trajectory of center of mass 



f decos(— f + e + 0) exp [— a(l — cose) 

< COS 9 > = : ; 

J aeexp[— a(l — cose)] 

. , f de cose exp [a cose] 

= sin^ p- f T ± 

J de exp [a cose] 

d f 

= sin</)— log{ / de exp [a cose]}. (51) 
da J 

By using the integration form of 0th modified Bessel 
function 



+ 1T 

Io(a) = — J de exp [a cose], 



(52) 



we obtain the mean velocity 



MU lo(a) 



(53) 
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while Iq/Iq is the odd function of a and so the signa- 
ture of B does not change the form of function. 

Note that from the asymptotic expansion of modified 
Bessel function, we obtain 



lim -pH±l. 

a— >±oo lo(a) 

Therefore we have 

<*>-> { i-£j*l}tf, 



(54) 



(55) 



in zero temperature limit, which is consistent with pre- 
vious result gll) , flU . 

The thermal expectation value of kinetic energy of cen- 
ter of mass is given as 



< Ekin > 



M 



<v 2 x + v 2 y > 



E + Muj c IU{1 



Where 



Eq 



M 



{U-u c lf 



}. (56) 



(57) 



essentially the same as the periodic potential energy in 
E e f f as shown in ([30| and (|31jl . Therefore if the temper- 
ature increases and the energy of relative motion (pendu- 
lum) comes up, mean value of potential energy is always 
finite and so as the kinetic energy. 

This temperature dependence of kinetic energy is ab- 
sent when the magnetic field disappeared. Such a de- 
pendence shows the energy transition between center of 
mass motion, relative motion and heat bath as we have 
expected. If the temperature of heat bath decreases, the 
energy of relative motion transfers to the bath, and en- 
ergy of center of mass transfers as well. 




FIG. 6: < Ekin > —Eo in the unit of Mu c lU as a function of 
dimensionless temperature T/T* 



VI. SUMMARY 



The figure 6 shows < Ekin > —Eq in the unit of 
MtUclU as the function of dimensionless temperature 
T/T*, where 



B | P U 



This function is monotone increasing function from to 
1 for the temperature increases from to oo. Therefore 
the mean kinetic energy has the asymptotic values as 



< Ehin > 



E 



M 



(U - u c lf for T ^ 0, (58) 



_> M-U^ + E^jf forT^oo. (59) 

The reason why the kinetic energy is finite under the 
infinite temperature is as follows. The kinetic energy is 



Here we have shown the new mechanism that the mag- 
netic field connects two dynamical systems, and if one 
system has the dissipation, another system becomes also 
dissipative only when the magnetic field is applied. This 
mechanism might also work in many other physical sys- 
tems under the magnetic field. Another interesting point 
is that when the electric dipole is placed in magnetic field, 
we have a constant with dimension of velocity. Then 
some kind of "resonance" phenomenon occurs and we 
have critical value of magnetic field which makes lowest 
kinetic energy vanish. 



B* 



MU 



(60) 



where U is a value of center of mass velocity when the 
magnetic field is absent. 
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